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SUMMARY 


PURPOSE 

The  purpose  of  this  program  was  to  develop  an  analytical  model 
describing  the  atmospheric  absorption  of  laser  energy  and  the  subsequent 
kinetic  processes  of  energy  exchange.  The  final  result  ol  this  effort  is 
the  quantitative  determination  of  the  atmospheric  temperature  change 
resulting  from  this  absorption. 

The  major  effort  has  been  concentrated  on  the  DF  laser  system. 

Eieht  representative  strong  laser  transitions  were  selected.  The  atmos¬ 
pheric  absorbers,  energy  levels  and  absorption  cross  sections  have  been 
specified.  The  analytical  kinetic  model  has  been  developed  and  rate 
constants  have  been  selected  for  use  in  the  calculation.  With  a  simple 
square  wave  pulse  shape,  the  model  has  been  used  to  give  an  assessment 
of  the  kinetic  effects  of  the  absorbers  in  heating  the  atmosphere  both  at 
sea  level  and  at  an  altitude  of  12  kilometers.  Two  limiting  concentrations 
of  water  have  been  used  and  the  laser  energy  has  been  varied  over 
several  orders  of  magnitude.  The  results  are  summarized  below  and 
detailed  in  subsequent  sections.  The  sea  level  atmospheric  temperature 
change  resulting  from  absorption  of  CO  laser  energy  by  atmospheric 
water  is  also  given. 

CONCLUSIONS 

A  model  has  been  developed  to  describe  the  kinetics  of  absorption 
of  DF  and  CO  laser  radiation  by  the  atmosphere. 

DF  laser  radiation  is  absorbed  by  N2  (collision  induced  continuum 
band),  C02  (v3),  NzO  (2vL  and  V]_  +  2v2),  CH4  (2v4  and  v2  +  v4),  HDO  (vj) 
and  H20  (v2  and  2v2)  and  the  H20  continuum.  Most  of  the  absorption 
bands  from  the  ground  vibrational  level  are  well  documented  in  the 
literature,  but  considerable  discrepancies  exist  concerning  absorption  by 
H20  and  HDO.  No  information  is  available  regarding  absorption  from  the 
upper  vibrational  levels,  and  these  have  been  ignored.  Under  conditions 
when  the  absorbing  transition  is  bleached,  this  may  become  a  severe 
limitation  for  the  model.  CO  laser  radiation  is  absorbed  by  C02  and  H?0. 
Of  these  two,  H20  (v2)  is  overwhelmingly  the  most  important  band,  and 
this  band  is  well  documented  in  the  literature. 

The  kinetics  model  describes  the  vibrational  relaxation  of  these 
absorbing  species  in  considerable  detail.  Rotational  relaxation  is  shown 
to  be  unimportant  for  intensities  below  10  W/cm  .  Considerable 
information  exists  regarding  the  relaxation  of  N2,  C02  (v3),  CH4  (v2  and 
v4)  and  H20  (v2).  However,  very  little  useful  information  exists  regarding 
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the  relaxation  of  N2O  (ej  and  i^)  and  HDO  (vj).  Since  the  relaxation  of 
1IDO  is  presumed  (by  analogy  to  l^O)  to  he  very  fast,  these  rates  are  not 
too  critical.  However,  N2O  is  an  important  absorber  and  the  effects  of 
vibrational  relaxation  in  N2O  are  significant.  Thus,  the  order-of-magnitude 
uncertainties  which  exist  in  some  of  the  important  relaxation  processes 

represent  a  serious  limitation  in  the  model. 

Preliminary  results  have  been  obtained  for  the  absorption  of  CO  and 
DF  radiation  under  a  variety  of  conditions.  The  results  indicate  that  N2O 
will  bleach  at  intensities  of  10^  to  10^  W/cm^,  depending  on  the  pressure 
and  humidity.  Under  conditions  of  very  iow  humidity,  absorption  by  CH4 
will  cause  cooling  of  the  atmosphere  for  short  times.  Absorption  of  CO 
laser  radiation  by  H2O  (v2)  and  subsequent  heating  of  the  atmosphere  is 
much  stronger  than  any  of  the  — ocesses  affecting  DF.  Under  conditions 
of  low  humidity,  however,  V  -V  transfer  from  H2O  (v 2)  to  O2  will  slow  the 
rate  of  heating  of  the  atmosphere. 
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L  SPECTROSCOPY 


A.  SELECTION  OF  DF  LASER  TRANSITIONS 

In  a  DF  chemic?.l  laser  the  DF  molecules  are  formed  by  chemical 
reaction  in  high  vibrational  and  rotational  states.  As  they  lase  and  relax 
by  collisional  processes  they  cascade  down  to  lower,  thermally  populated 
levels.  Thus,  in  pulsed  DF  lasers  the  various  transitions  do  not  all 
lase  simultaneously.  Rather,  the  higher  vibrational  levels  tend  to  lase 
first  with  lower  levels  appearing  later  in  the  pulse.  Even  in  continuous 
wave  (CW)  lasers  the  relative  strength  of  the  various  lines  is  a  complex 
function  of  the  chemical  and  relaxation  processes  occuring  in  the  laser. 

The  various  DF  lines  are  absorbed  by  different  species  with  widely  vary¬ 
ing  absorption  cross  sections.  Therefore,  it  may  be  found  desirable  to 
suppress  certan  transitions,  even  though  they  may  be  strong  laser  lines, 
if  they  are  strongly  absorbed  and  adversely  affect  the  propagation  of  other 
lines.  However,  this  may  be  difficult  to  accomplish  in  DF  lasers  due  to 
their  intrinsic  high  gain  and  the  complex  nature  of  the  processes  control¬ 
ling  the  relative  strength  of  the  various  lines. 

The  eight  DF  laser  transitions  shown  in  Table  I  were  chosen  from 
the  experimental  work  on  pulsed  DF  lasers  of  Deutsch^ )  and  Basov, 
et  al.  V2)  Reference  1,  which  contains  tire  most  accurate  experimental 
measurements  of  line  position  (±  0.  08  cm'*),  was  used  to  determine  the 
transition  frequency,  and  the  relative  energy  data  of  Ref.  2  were  used  to 
select  lines  of  strong  output  in  pulsed  lasers.  Care  was  employed  to 
ensure  that  all  known  atmospheric  absorbers  were  included  for  these  lines. 

TABLE  I 


PULSED  DF  LASER  TRANSITIONS 


Laser  Line 
Identification 

La.’^r  Frequency 
(cm'  1 ) 

Relative 

Energy* 

2P  (8) 

2631.  09 

2.  37 

2P  (9) 

2605.  87 

2.  68 

2P  (10) 

2580.  16 

4.  24 

2P  (11) 

2553.  97 

3.  46 

3 P  (8) 

254b.  3? 

2.  74 

3 P  (9) 

2521.  81 

5.  38 

3P  (10) 

2496.  61 

3.  35 

3P  (11) 

2471.  34 

2.  26 

$ 

Energy  of  IP  (10)  transition  is  taken  as  1. 


B.  SELECTION  OF  ATMOSPHERIC  ABSORBERS  AND  ABSORPTION 
COEFFICIENTS  FOR  DF  LAE..R  LINES 

The  atmospheric  absorbers  of  importance  at  DF  laser  wavelengths 
are  N2O,  CH4,  H?0,  HDO  and  CO2.  This  is  confirmed  by  both  absorption 
cell  experiments^)  and  theoretical  calculations  ba^ed  on  experimental 
data.  (4)  Reference  4  is  an  Air  Force  Cambridge  Research  Laboratory 
(AFCRL)  compilation  of  the  molecular  spectroscopic  parameters  for  a 
number  of  infrared-active  molecules  (including  those  listed  above)  which 
occur  naturally  in  the  ter restrial  atmosphere.  Parameters  included  in 
the  compilation  for  each  absorption  line  are:  frequency,  intensity,  half¬ 
width,  energy  of  the  lower  state  of  the  transition,  vibration  and  rotational 
identifications  of  the  upper  and  lower  energy  states,  an  isotopic  identifica¬ 
tion,  and  a  molecule  identification.  We  have  written  a  code  which  sorts 
the  AFCRL  data  according  to  absorbing  molecules  and  vibrational  transition 
over  a  given  frequency  interval  and  calculates  the  absorption  cross  section 
for  the  absorbing  molecule  at  the  laser  frequency  of  interest. 

The  sea  level  absorption  cross  sections  used  for  the  atmospheric 
absorbers  at  DF  laser  frequencies  are  given  in  Table  II  The  absorption 
by  N2  on  the  collision  induced  vibrational  continuum  band  centered  around 
4.  3  ji  is  not  included  here,  but  is  a  well  documented  phenomenon.  Calcu¬ 
lations  indicate  that  the  N2  continuum  is  a  strong  absorber  at  the  longest 
wavelengths  at  sea  level.  Following  McClatchyf^)  the  effect  of  absorption 
due  to  the  water  continuum  was  not  included.  This  effect  will  be  investi¬ 
gated  in  a  later  phase  of  this  effort.  At  higher  altitudes  the  cross  sections 
are  reduced  according  to  the  appropriate  pressure  dependence.  The 
temperature  dependence  is  not  included  at  this  time.  (See  Section  III, 
RESULTS.  ) 

The  N2O  cross  sections  computed  are  in  excellent  agreement  with 
the  absorption  cell  measurements  of  Ref.  3  for  all  the  laser  transitions 
considered.  For  CH4  the  agreement  between  these  cross  sections  and 
the  absorption  cell  measurements  of  Ref.  3  are  not  as  good.  In  all  cases 
we  have  used  the  cross  sections  computed  from  the  AFCRL  data,  which 
average  about  a  factor  of  3  to  4  lower  than  the  absorption  cell  measure¬ 
ments.  The  effects  of  any  errors  thus  introduced  on  the  calculated  atmos¬ 
pheric  temperature  rise  are  small,  because  in  no  case  is  CH4  the  major 
absorber.  For  CO2,  where  the  absorption  is  due  to  lines  of  the  4.  3  p.  band 
(hundreds  of  wave  numbers  removed  from  the  laser  frequency)  we  have  used 
the  experimental  results  of  Ref.  3.  These  fall  between  the  values  calculated 
form  the  AFCRL  data  using  a  Lorentz  line  shape  (higher  than  Ref.  3)  and 
the  best  experimental  line  shape  fit  to  CO2  absorption  data(b)  (lower  than 
Ref.  3). 

The  absorption  cross  section  is  calculated  from  the  integrated 
absorption  coefficient  assuming  a  Lorentz  line  shape  using  the  relation: 

S  _ Av _  n» 

CT  “  7T  ,  .2  2 

(v-  v  )  +  Av 

'  o 
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TABLE  II 

SEA  LEVEL  ABSORPTION  CROSS  SECTIONS  AT 
DF  LASER  WAVELENGTHS 
{ cm^/  molecule) 


DF  Laser 
Line 


2P  (11) 
2P  (10) 
2P  (9) 
2P  (8) 
3P  (11) 
3  P  (10) 
3  P  (9) 
3P  (8) 


N2° 

HDO 

H2° 

ch4 

1.0-20 

-  23 

3.  7  x  10 

8.  0  x  10“29 

- 

2..9-20 

-22 

1.  77 

1. 4  x  10"  28 

- 

3.  57“22 

1. 53“21 

1.  0  x  10"27 

-23 

1.5 

5.  O’  22 

6.  8'26 

2.  0"22 

4.2-21 

6.  7  x  10-23 

i.r25 

- 

2.  39“22 

3  x  1 0*  24 

9.3-26 

- 

4.  7- 22 

-23 

1.2x10 

8.  0  x  10"28 

- 

2.  9- 20 

-23 

5.  7  x  10  " 

7.  4  x  10’28 

2.5-22 

3.  15 


where  a  is  the  cross  section  in  cm^/pa  rticle,  S  (the  integrated  absorption 
cross  section  cm- 1 /particle- cm_2)  =  j  a  (v  )  dv  ,  Av  is  the  half-width  in 
cm"  1  and  (v  -v>  Q)  is  the  difference  in  wave  number  between  the  laser  line 
frequency  and  the  center  of  the  absorbing  line.  This  procedure  gives  too 
high  a  value  for  CO2,  since  CO?  is  known  to  be  sub- Lorent/ian  far  from 
line  center.  Reference  6  reports  measurements  of  CO2  absorption  both 
in  pure  CO2  and  in  mixtures  of  CO2  with  N2  and  O2.  The  absorption  due 
to  the  wings  of  the  strong  lines  in  the  v  3  band  centered  at  2349  cm- 1  is 
much  less  than  that  calculated  with  the  Lorent/  shape.  The  data  in  N2  and 
O2  are  fit(6)  with  the  following  equation,  which  retains  the  Lorent/.  pressure 
dependence  but  requires  a  nearly  exponential  modification  of  the  frequency 
dependence. 


S  Av 

a  -  —  - = - 

7T  .  ,2  ,  . 

(v  -v  )  +  Av 

'  o' 

Use  of  this  relation  gives  values  for  CO£  absorption  cross  sections  at  the 
DF  laser  frequencies  of  interest  which  are  lower  than  those  of  Ref.  3.  The 
cross  sections  for  CO2  given  in  Table  II  are  consistent  with  the  measure¬ 
ments  in  Ref.  3.  It  should  be  noted  that  any  errors  introduced  by  uncer¬ 
tainties  in  the  value  of  the  CO2  cross  sections  would  be  significant  only  at 
high  altitudes  where  CO2  becomes  a  significant  absorber  (see  Fig.  14). 

At  sea  level  CO2  only  a  very  minor  contributor  to  the  atmospheric 
temperature  rise. 

For  atmospheric  water  vapor  the  situation  is  complicated  by  the 
fact  that  both  H2O  and  HDO  are  absorbers,  with  HDO  such  a  strong  absorber 
that  it  is  significant  even  in  its  natural  isotopic  abundance  of  300  ppm  of 
H2O.  The  calculated  cross  sections  for  H2O  and  HDO  obtained  from  the 
AFCRL  data  and  the  experimental  results  of  Ref.  3  aie  in  sharp  disagree¬ 
ment  for  individual  laser  lines.  Because  the  experimental  absorption  cross 
sections  were  obtained  by  subtracting  the  very  large  absorption  due  to  D?0 
from  the  measured  absorption,  we  have  rather  arbitrarily  assumed  the 
error  to  be  in  the  experiment  and  have  chosen  to  use  the  values  calculated 
from  the  AFCRL  data  which  for  H2O  are  in  agreement  with  data  of  Bates.  (?) 
It  turns  out  that,  when  the  laser  power  is  evenly  distributed  over  all  the 
laser  lines  with  equal  powers,  either  set  of  cross  section  values  gives  the 
same  temperature  rise.  This  is  because  the  sum  over  all  lines  of  the 
absorption  cross  sections  for  HDO-H2O  happens  to  be  the  same  for  both 
the  AFCRL  and  Ref.  3  cases.  It  should  be  noted,  however,  that  for  indi¬ 
vidual  laser  lines  the  values  are  very  different  and  that  any  calculations 
of  the  temperature  rise  due  to  H2O-HDO,  which  use  less  than  all  of  the 
eight  laser  lines  of  Table  I  or  employ  different  power  levels  on  the  lines, 
may  be  sensitive  to  which  set  of  cross  section  values  is  chosen. 


-0.  46  (1  v  -v  -  5) 
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C.  SELECTION  OF  CO  LASER  LINES 

In  contrast  to  DF  chemica]  lasers,  the  CO  molecules  in  electric 
CO  lasers  are  first  excited  to  relatively  low  vibrational  levels  by  the 
discharge.  I  hey  are  then  pumped  to  rather  high  vibrational  levels  by 
mtermolcular  V-V  transfer  processes.  Due  to  anharmonic  effects  the 
vibrational  distribution  function  is  generally  strongly  distorted  from  a 
Boltzmann  distribution  in  a  manner  which  favors  the  higher  vibrational 
levels.  As  a  result,  the  highest  gain  is  usually  observed  in  relatively 
high  levels  (typically  !  <  v  £  15).  Recently,  lasing  has  been  observed 
on  the  lower  levels  (down  to  the  1-0  transition)  in  high  gain-low  loss 
systems. 

As  we  shall  see,  the  higher  vibrational  levels  which  emit  at 
longer  wavelengths,  are  more  strongly  absorbed  by  tlu  H20  (v?)  6.  3-n 
band.  In  addition,  the  kinetics  of  the  lower  levels  take  place  more  rapidly 
offering  the  possibility  of  greater  specific  power.  For  these  reasons  it 
may  oe  desirable  to  suppress  lasing  on  the  high  vibrational  levels  (with  a 
grating  or  passive  H20  absorber,  for  example)  and  force  lasing  on  the  low 
vibrational  levels.  Experiments  to  demonstrate  the  feasibility  of  this  are 
currently  underway  at  AERL. 

To  assess  the  relative  advantages  of  low  and  high  vibrational  levels 
.Vs  "e"eSSary  t0  conslder  a  wil  e  range  of  transitions.  The  positions  of 
the  CO  laser  lines  used  in  calculating  the  results  presented  here  are  given 
m  Table  III  together  with  absorption  coefficients  and  cross  sections  for 
absorption  by  H2Of  the  only  important  absorber.  The  wavenumber  positions 
or  CO,  listed  in  Table  III,  were  calculated  from  spectroscopic  constants 
given  by  Mantz.  (oa,  b)  These  should  be  good  to  ±0.  01  cm-1.  The  wave- 
numbers  agree  well  with  those  used  by  Long(9)  and  McClatchy.  (10) 


D. 


SELECTION  OF  CO  ATMOSPHERIC  ABSORBERS  AND  ABSORPTION 
COEFFICIENTS 


Survey  calculations  were  made  of  absorption  of  CO  laser  lines  by 

at2S’  Nz°  a?d.  CH4-  Compared  with  H20,  the  absorption  of  O5, 

%Z^.and  CH4  was  negligible.  The  absorption  of  the  6P  (9)  CO  line  at 
cm_",  hy  C°2  was  als°  negligible.  For  the  2P  (9)  CO  line  at 
2UH1.  26  cm”  ,  the  C02  absorption  was  down  from  that  due  to  H?0  by  a 
factor  of  6.  (The  sub-Lorentz  C02  line  contour  of  Winters,  et  al.  ,(°)was 
used.  )  In  the  remainder  of  this  work  we  concern  ourselves  only  with 
absorption  by  H20. 

Fxtensive  calculations  of  absorption  of  CO  laser  radiation  by 
atmospheric  H^O  have  been  made  by  McClatchey.  U0,  H)  However,  this 
work  has  two  limitations.  First,  the  Lorentz  contour,  Eq.  (1)  was  used 
for  the  HoO  spectral  absorption  coefficient.  The  experiments  of  Long, 
et  al.  ,  )  show  that  the  absorption  of  CO  laser  lines  in  the  wings  of  ILO 

lines  is  typically  40%  greater  than  given  by  Eq.  (1).  To  fit  their  data  they 
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assumed  that  the  absorption  coefficient  follows  the  Lorent/  contour  for 
|  v  -v’ol  <  wm,  but  that  outside  this  range  it  is  given  by 


ct(v  )  =  - 


(v  r  av 

m 


7!  r  .2  .  2")  rn  1  0 

pm’  +Al'Jr-vo|  • 


Equations  (1)  and  (3)  agree  a'.  |v-v0|  =  vm.  The  best  fit  was  obtained  using 
v’m  =  3  Av  and  m  =  1.  77. 


McClatchey  and  Long  considered  absorption  by  H^O  lines  whose 
centers  lie  within  20  or  25  cm- 1  ,  respectively,  of  the  CO  laser  line.  This 
amounts  to  a  cut-off:  a  (v)  =  0,  v  >  vmaj?  or  v  <  vmir|.  However,  we  find 
by  allowing  ^min  =  vo  -  350  cm  “  that  in  some  cases  half  the  absorption  is 

due  to  lines  150  cm-1  distant.  The  true  contour  in  the  far  wings  of  H2O 
lines  cannot  be  determined  from  the  available  data.  Unfortunately,  H^O 
does  not  form  sharp  band  heads  as  does  CO2,  whose  far  wing  line  profile 
may  be  determined  more  reliably.  We  consider  extrapolation  of  Eq.  (3) 
less  arbitrary  than  cutting  it  off.  Using  the  extended  line  profile  we  find 
that  Long'  s  data  may  be  fitted  using  vm  =  3  L  and  m  ~  1.  88.  Values  o  m 
for  various  individual  lines  ranged  from  1.  85  to  1.91.) 


In  mixtures  of  water  vapor  and  air  the  H2Q  half-width  is  given  by 


Av  =  %  Pt  U  +(B-l)xl  (4) 

where  aQ  is  the  half-width  per  atmosphere  of  air,  pt  is  the  total  pressure 
in  atmospheres,  X  is  the  mole  fraction  of  H?0  in  the  mixture  and  B  is  the 
self-broadening  coefficient.  Long,  et  al.  ,  (9)  find  B  =  8  to  27.  We  used 
the  value  B  =  5  from  Ref.  12  The  difference  is  insignificant  because 
X  <  0.  01  in  cases  of  interest. 

Our  results  are  presented  in  Table  IIL  Equation  (3)  was  used  with 
vm  =  3Av  and  m  =  1.  88.  The  cut-offs  were  ‘'min  -  vo  ~  350  cm*l  and 
’'max  =  +  50  cm'*,  {Absorption  at  larger  wavenumbers,  farther  away 

from  the  center  of  the  6.  3-p.  H2O  band  is  negligible.  )  The  absorption 
coefficients  are  larger  than  those  scaled  from  McClatchey'  s  calculations(  1  0) 
for  PH2O  =  3.  3  torr  (mid- latitude  winter  case),  due  to  the  use  of  the 

"super-Lorentz"  profile  and  its  extension  farther  into  the  wings  of  the  H?0 
lines.  There  are  no  Ohio  Stated)  experimental  or  theoretical  results  for 
many  of  these  lines.  Their  work  has  generally  been  on  lines  at  smaller 
wavenumbers  (higher  CO  vibrational  levels)  located  nearer  the  center  of 
the  6.  3-fj.  baud  of  H2O,  for  which  the  attenuation  coefficients  are  typically 
larger. 

Table  III  shows  the  cross  sections  we  have  used  for  water  based 
on  the  above  discussion.  For  comparison  we  have  shown  the  differences 
between  our  results  and  those  of  McClatchey^®)  in  terms  of  absorption 
coefficients.  (Note  added  in  proof:  A  recent  work^  ^)  reports  absorption 
coefficients  in  good  agreement  with  the  present  ones,  ) 


11  KINETICS  OF  ABSORPTION  OF 
D,  AND  CO  LASER  RADIATION 


The  absorption  of  DF  laser  radiation  is  complicated  by  the  number 
of  atmospheric  species  which  absorb  this  radiation  and  by  the  fact  that  the 
upper  level  of  the  absorbing  transition  is  frequently  a  high-lying  vibrational 
level,  as  in  N2O  and  CH^.  Thus,  to  describe  the  absorption  kinetics  it  is 
necessary  to  include  a  large  number  of  levels.  A  complete  diagram  of 
all  the  important  vibrational  levels  of  the  atmospheric  species  important 
for  DF  and  CO  absorption  kinetics  is  shown  in  Fig.  1.  It  is  important  to 
keep  in  mind  that  each  vibrational  level  shown  actually  consists  of  a 
manifold  of  vibration-rotation  states. 

A.  Rotational  Relaxation 

Because  of  the  large  number  of  levels,  it  is  necessary  to  make 
approximations  in  which  groups  of  levels  are  lumped  together  and  treated 
as  a  single  ievel.  To  begin  with,  it  will  be  assumed  that  all  the  rotational 
levels  belonging  to  a  given  vibrational  level  remain  in  equilibrium  with 
one  another  at  the  translational  temperature  T.  This  assumption  will  be 
valid  if  the  rotational  relaxation  time  Tr  is  short  compared  with  the  time 
for  a  radiatively  induced  transition,  that  is,  if 

$  a  7  D  «  l ,  (9 ) 

v  v  R 


where  is  the  flux  of  laser  photons,  and  av  is  the  absorption  cross  section. 
Both  the  rotational  relaxation  time  and  the  absorption  cross  section  depend 
upon  the  pressure  throughthe  molecular  collision  frequency  vc.  Thus,  the 
adequacy  of  assumption  (5)  will  depend  on  the  pressure.  For  a  purely 
Lorentz  broadened  line,  the  cross  section  at  the  frequency  v  is  given  by 
Eq.  (1).  The  line  half -width  Av  is  proportional  to  the  pressure  and  may 
be  represented  by  the  relation 


Av  =  ap, 


where  p  is  the  pressure. 

The  rotational  relaxation  time  is  more  difficult  to  quantify  precisely. 
Since  rotational  relaxation  involves  a  great  many  levels,  it  is  inevitably  a 
complex  process.  For  example,  an  extremely  short  laser  pulse  (shorter 
than  the  time  between  molecular  collisions)  will  affect  only  molecules  in 
specific  quantum  states.  The  population  of  the  perturbed  states  will  then 


Preceding  page  blank 
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relax  by  reorientation  of  the  remaining  molecules  having  the  same  total 
angular  momentum  and  by  transitions  from  nearby  rotational  energy  levels. 
The  appropriate  relaxation  time  for  this  case  corresponds  to  the  gain 
relaxation  time  measured  following  a  Q-switched  laser  pulse.  (^)  Longer 
laser  pulses  will  pertrub  the  populations  of  all  the  rotational  levels  around 
the  level  from  which  the  transition  originates;  in  this  case  the  time  for  the 
system  to  relax  to  equilibrium  will  be  considerably  longer.  Moreover, 
experimental  evidence  indicates  that  the  thermally  averaged  relaxation 
time  (measured  acoustically)  may  correspond  to  as  small  as  a  few  collisions 
or  as  large  as  30  collisions  or  more  (for  rotaticnally  light  molecules,  such 
as  HC1,  at  elevated  temperatures).  (15)  Moreover,  it  has  been  found 
both  theoretically  and  experimentally  that  for  molecules  in  high  rotational 
levels  the  rotational  relaxation  time  increases  exponentially  with  rotational 
energy.  (^) 

To  examine  the  importance  of  rotational  relaxation  it  is  convenient 
to  express  the  relaxation  time  in  the  form 


PTr  =  P. 


(6) 


where  p  is  independent  of  pressure.  Equation  (4)  may  then  be  expressed 


a  3 


$  a  r  n  =  4> 
v  v  R  V  7T 


(v  -  vo)2  +  (ap)2 


«  1.  (7) 


When  the  laser  line  lies  in  the  wing  of  the  absorber  (  v  -  v  |/ap  »  1) 
we  obtain  the  inequality 


<I>  « 

v 


ir  (v  -  v  ) 
o 

Q  P  S 


(8) 


The  right-hand  side  is  a  function  only  of  the  laser  frequency  and  the  absorbing 
level,  and  is  independent  of  the  pressure.  The  limiting  fluxes  for  the  most 
important  DF  laser  lines,  calculated  using  this  expression,  are  summarized 
in  Table  IV.  At  sufficiently  high  pressures  (a  p/ |  v  -  v0  |  »  1)  the  laser  line 
will  lie  within  the  Lorentz  width  of  the  absorber.  The  limiting  flux  $v  will 
then  increase  as  p2.  At  sufficiently  low  pressures,  of  course,  the  line  will 
become  Doppler  broadened  and  the  cross  section  a ^  will  be  independent  of 
pressure.  The  limiting  flux  $v  at  which  rotational  relaxation  becomes 
important  will  then  be  proportional  to  the  pressure;  at  a  sufficiently  low 
pressure  rotational  relaxation  will  always  become  important.  However, 
as  shown  in  Table  IV,  this  occurs  only  at  very  low  pressures  (<,  0.  01  atm). 

When  rotational  relaxation  becomes  important  the  atmosphere  will 
begin  to  "bleach",  that  is,  molecules  will  be  removed  from  the  absorbing 
level  and  deposited  in  the  upper  level  faster  than  they  can  be  relaxed  by 
collisions  and  the  atmosphere  will  become  transparent.  Thus,  by  ignoring 
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TABLE  IV 


DF  LASER  INTENSITY  AT  WHICH  ROTATIONAL 
RELAXATION  BECOMES  IMPORTANT 


DF  Laser 
Transition 

Abso  -bing 
Species 

Relaxation 

Time 

(collisions ) 

Lase  r 
Intensity 
(W  /cm^) 

Pressure 

P.ange 

(atm) 

2P  (8) 

CH4 

10 

8  x  109 

0.  01  <  p  <  0.  2 

2P  (10) 

n2o 

2 

4  x  109 

0  <  p  <  1 

2P  (11) 

n2o 

2 

IxlO10 

0  <  p  <  3 

3P  (8) 

IS  O 
c 

2 

1  x  109 

0.  02  <  p  <  0.  04 

Notes: 


At  lower  pressures  the  absorbing  lines  become  Doppler  broadened, 
and  the  intensity  at  which  rotational  relaxation  becomes  important  decreases 
linearly  with  the  pressure.  At  higher  pressures  the  laser  line  falls  within 
the  uorentz  width  of  the  absorbing  line  and  the  intensity  at  which  rotational 
relaxation  becomes  important  increases  as  the  square  of  the  pressure. 

Laser  transitions  and  absorbing  species  not  indicated  in  Table  IV 
show  no  rotational  relaxation  effects  at  laser  intensities  below  1010  W/cm2. 
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rotational  relaxation  we  place  an  upper  bound  on  the  extent  of  absorption. 

The  results  shown  in  Table  IV  indicate  that  rotational  bleaching  should  not 
be  important  for  DF  laser  intensity  levels  below  10  ®  W/cm^.  Although 
a  few  isolated  lines  may  be  bleached  at  this  power  >evel,  their  contribution 
to  the  total  absorption  on  the  vibrational  bands  to  which  they  belong  will 
not  be  significant.  It  should  be  emphasized,  however,  that  this  result 
cannot  be  simply  extrapolated  to  all  IR  lasers.  All  the  important  absorbers 
for  DF  radiation  have  rather  small  absorption  cross  sections  because  either 
the  absorbing  bands  correspond  to  vibrational  overtones,  (as  in  the  case 
of  N2O  and  CH4),  or  the  laser  line  lies  in  the  far  wings  of  the  absorbing 
line  (as  in  CO2),  or  the  rotational  structure  of  the  band  reduces  the  contribu¬ 
tion  of  individual  lines  (as  in  HDO  and  H2O)  For  laser  lines  overlapped 
by  strong  absorbing  lines,  rotational  relaxation  may  become  important. 


B.  VIBRATIONAL  RELAXATION 


Even  if  rotational  equilibrium  is  assumed,  there  remain  a  great 
many  vibrational  levels  of  importance  for  the  absorption  kinetics,  as  indica¬ 
ted  in  Fig.  1.  Although  it  is  not  difficult  to  handle  this  number  of  levels  on 
the  computer,  it  is  appropriate  to  make  some  further  simplifications  since 
liiere  do  not  exist  enough  kinetic  rate  data  to  enable  one  to  describe  the 
system  in  complete  detail.  The  kinetics  of  the  important  atmospheric 
absorbers  are  discussed  in  the  following  sections. 

1.  Kinetics  of  Absorption  of  DF  Laser  Radiation  by  N2 

DF  laser  radiation  is  absorbed  by  N2  on  the  collision  induced 
vibrational  continuum  band  centered  around  4.  3  fi: 


N2  (  v  =  0)  +  hv  +  M  -  N2  (  v  =  1)  +  M  +  AE.  (9) 


The  excess  energy  AE  =  hv  -  Evib  ~  250  cm-1  appears  immediately  as 
translational  and  rotational  heating  of  the  N2  and  M  molecules. 

The  relaxation  of  N2  (v  =  1 )  by  the  dominant  atmospheric  species 
(N2,  O2,  H2O)  has  been  extensively  studied  as  a  result  of  interest  inthe 
CO2  -  N2  laser.  T-V  deactivation  of  N2  (v  =  I )  by  N2  and  O2  and 

V-V  transfer  to  O2  are  extremely  slow.  The  most  important  relaxa¬ 
tion  path  is  provided  by  H2O.  It  has  not  been  established  directly  *hether 
the  quenching  of  N2  by  H2O  proceeds  by  a  T-V  (or  vibration  to  rotation) 
process,  or  by  V-V  transfer  to  the  V2  mode  of  H2O, 


N2  (v  =  1)  +  HzO  (000)  -  N2  (  v  =  C)  +  HzO  (010)  +  735  cm'1  (10) 


However,  the  subsequent  relaxation  of  H2O  (010)  is  so  fast  that  for  practical 
purposes  it  does  not  matter. 


These  processes  are  summarized  in  Fig.  2.  The  time  scales  for 
the  various  processes  in  this  and  the  following  sections  are  based  on  an 
atmosphere  composed  cf 


(2  x  10^  molecules/cm^), 
(5  x  10^  molecules/cm^), 
f^O  (3  x  1017  molecules/cm^), 


at  a  temperature  of  300°K.  The  amount  of  energy  which  appears  as 
thermal  (translational  +  rotational)  energy  at  each  step  is  also  indicated. 

We  see  that  under  the  above  conditions,  absorption  on  the  N2  continuum 
band  results  in  a  small  amount  cf  essentially  instantaneous  heating  followed 
by  slow  heating  on  the  N2  relaxation  time  scale,  which  is  on  the  order  of 
8  x  10-4S. 

2.  Kinetics  of  Absorption  of  DF  Laser  Radiation  by  CO2 

DF  laser  radiation  is  absorbed  by  CO2  on  the  V3  fundamental  band 
centered  at  4.  3  /J  . 


C02  (000)  4  hv  -  COz  (001)  +  AE. 


(11) 


The  absorption  takes  place  principally  in  the  far  wings  of  the  lines  corres¬ 
ponding  to  low,  thermally  populated,  rotational  levels.  Thus  most  of  the 
excess  energy  AE  (-'-250  cm'*)  appears  instantaneously  as  translational 
energy  of  the  CO2  and  collisional  broadening  molecules. 

The  relaxation  of  CC>2  (001)  by  the  dominant  atmospl  ric  species 
(N?,  O2,  H2O)  has  been  extensively  studied  as  a  result  of  interest  in  the 
CCTj  -  N2  laser.  (17)  Relaxation  takes  place  principally  via  intermolecular 
V-V  transfer  to  N2. 

C02  (001)  +  N2  (v  =  0)  -  C02  (000)  +  N2  (v  =  1)  -  18  cm'*  (12) 


The  subsequent  relaxation  of  N2  (  v  =  1)  is  described  above. 

These  processes  are  summarized  in  Fig.  3,  where  the  relaxation 
times  correspond  to  the  same  atmospheric  conditions  as  in  Fig.  2.  Aside 
from  the  small  amount  of  heat  released  by  the  absorption  process,  the 
overwhelming  effect  of  absorption  by  CO2  is  to  heat  the  gas  on  a  slow 
time  scale  (~10"^s)  corresponding  to  the  relaxation  of  the  N2  (v  =  1)  level. 
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E  ~  250  cm1  (TYPICAL) 
^ -  y  =1 


n  of  DF  Laser  Radiation  by  N£ 
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Kinetics  of  Absorption  of  CO  and  DF  Laser  Radiation  by  HpO 

CO  laser  radiation  is  absorbed  by  H2O  on  the  v  2  fundamental 
transition, 


HzO  (000)  +  hv  -  HzO  (010)  +  AE.  (13) 

The  excess  energy  (AE  ~  400  cm'l)  appears  as  translational  and  rotational 
energy  of  the  H2O  molecule,  and  is  assumed  to  be  thermalized  instantly. 

the  absorption  of  DF  radiation  by  H2O,  while  somewhat  controver¬ 
sial,  as  mentioned  above,  is  alleged  to  occur  on  the  same  vibrational 
transition.  The  relatively  large  excess  energy  (AE~900  cm'  for  DF 
lasers)  appears  as  translational  and  rotational  energy  of  the  H20  molecule, 
and  is  presumed  to  be  thermalized  instantaneously. 

The  relaxation  of  the  H2O  (v  2)  mode  has  been  studied  experimentally 
for  a  long  time,  and  while  considerable  uncertainty  remains  for  sorm  of 
the  important  rates,  a  broad  understanding  is  possible.  (17)  The  most 
important  processes  for  the  relaxation  of  H2O  (v->)  are  near  resonant 
V-V  transfer  to  O2, 


1^2^  (010)  +  C>2  (v  =  0)  -*  H2O  (000)  +  O2  (v  =  1)  +  38  cm‘\  (14) 
and  T-V  quenching, 


HzO  (010)  +  M  -*  H20  (000)  +  M  +  1595  cm'l,  (15) 

where  M  =  N2,  02,  and  (most  importantly)  H20  itself.  There  is  consider¬ 
able  scatter  in  the  data  on  the  T-V  quenching  of  H20  (v2)  by  H20,  but  the 
rate  is  apparently  quite  fast,  on  the  order  of  (2  x  10'12)  -  (5  x  1  0 “ 1 1 )  cm3/s, 
(about  10  to  100  collisions).  The  most  recent  data  favor  the  fastest  rate. 

The  data  on  V-V  transfer  to  O2  (reaction  (14))  are  somewhat  confused 
due  to  apparent  errors  in  most  of  the  papers.  The  correct  rate  seems 
to  be  about  1.  5  x  10"12  cm^/s.  The  energy  which  is  transferred  to  O2 
relaxes  relatively  slowly  by  T-V  quenching  of  O2  by  H20, 

C>2  (v  =  1)  +  H20  (000)  -*■  02  (v  =  0)  +  H20  (000)  +  1556  cm~l,  (16) 

and  by  transfer  back  to  H2O  followed  by  quenching  of  the  H20. 
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1  hese  processes  are  summarized  in  Fig.  4,  where  the  times  indicated 
correspond  to  the  atmospheric  conditions  of  Figs.  2  and  ^  Ihe  absorption 
of  a  laser  photon  leads  to  some  immediate  heating  (for  a  time  on  the  order  of 
the  rotational  relaxation  time)  as  indicated.  T-V  quenching  of  the  HpO 
leads  to  further  very  rapid  heating  (with  a  time  scale  on  the  order  of  10~7 
to  10'^s.  The  energy  that  is  transferred  to  O2  produces  heating  on  a 
somewhat  larger  time  scale,  on  the  order  of  lO'^s.  The  extent  of  heating 
occurring  on  the  fast  and  slow  t'me  scales  depends  on  the  relative  rates 
of  T-V  quenching  and  V-V  transfer  to  0£-  If  the  faster  HpO  O2  V -V 
transfer  rate  and  or  slower  H2.O/H2O  T-V  quenching  rate  is  correct,  or 
if  the  humidity  is  low  (as  it  will  sureiy  be  at  higher  altitudes)  then  the  bulk 
of  the  H2O  vibrational  energy  will  be  transferred  to  O^.  This  will  shift 
the  heating  to  significantly  longer  times.  These  uncertainties  become 
even  more  important  for  the  absorption  of  CO  laser  radiation,  since 
Il20(v2)  Is  dominant  absorber. 

4.  Kinetics  of  Absorption  of  DF  Laser  Radiation  by  N2O 

DF  laser  radiation  is  absorbed  by  N2O  principally  on  the  v2  overtone 
transition, 

NzO  (00°0)  +  hv  N20  (20°0)  4  AE  (17) 

and  the  combination  band  transition 

fi  f 

NzO  (00  0)  +  hv -»  NzO  (12  0)  +  AE  (18) 

A  small  amount  of  absorption  takes  place  on  the  corresponding  transitions 
from  the  N20  (01 1 0 )  level,  but  this  may  be  ignored. 

Most  of  the  available  data  on  N2O  vibrational  relaxation'  '  '  were 

obtained  in  pure  N2O  and  N2O/N2  mixtures,  and  do  not  shed  much  light 
on  the  most  important  processes  in  the  atmosphere.  To  identify  the  most 
important  processes  we  note  that  N2O  is  very  similar  to  CC^.  Both  are 
linear  triatomic  molecules  with  comparable  masses.  It  is  important  to 
note,  however,  that  mode  v^  is  IR  active  (that  is,  it  possesses  a  dipole 
moment)  in  N2O,  whereas  it  is  not  IR  active  in  C02-  This  is  because  CO2 
is  a  symmetric  molecule  whereas  N2O  is  not,  having  the  configuration 
N-N-O.  It  is  also  important  to  remember  that  the  levels  10°0  and  02°0 
are  strongly  coupled  by  Fermi  resonance  in  CO2  but  not  in  N2O  Available 
data  indicate  that  modes  vi  and  v2  of  C02  are  rapidly  coupled  by  collisions 
and  may  be  assumed  to  relax  together,  but  that  mode  V3  is  only  weakly 
coupled  to  the  other  two  and  relaxes  by  itself.  All  the  data  concern  mixtures 
rich  in  CO2,  so  that  modes  vj  and  V2  may  be  coupled  by  intermolecular 
V-V  transfer  processes  of  the  type 
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(19) 


CO  (10°0)  +  COz  (00°0)-*  2  COz  (01*0) 


or  -  C02  (00°0)  +  C02  (02C0). 


It  has  also  been  suggested(22-23)  that  in  dilute  mixtures  coupling  via 
intromolecular  processes, 


C02  (10°0)  +  f*l  -*  C02  (02^0)  +  M  (?0\ 

where  M  =  CO2,  is  rapid  due  to  the  Fermi  resonance  which  mixes  the 
levels  (10°0)  and  (02°0). 

There  is  no  corresponding  evidence  of  rapid  coupling  of  modes 
vj  and  v2  of  N20  in  the  atmosphere.  Intermolecular  V-V  processes  similar 
to  (19)  are  unimportant  due  to  the  high  dilution,  and  modes  vj  and  v?  are 
not  mixed  by  Fermi  resonance.  The  acoustical  data  have  been  interpreted 
in  terms  of  the  simultaneous  relaxation  of  modes  and  v2,  but  the  data 
are  not  compelling.  For  lack  of  evidence  we  shall  make  the  same  approxi¬ 
mation. 


The  evidence  does  show,  however,  that  mode  V3  is  only  weakly 
coupled  to  modes  vj  and  v2.  Results  obtained  in  N_20/N2  mixtures  indicate 
that  quenching  of  mode  V3  by  N2  requires  about  10  collisions,  which  is 
only  about  a  factor  of  two  faster  than  the  corresponding  process  in  COt. 

It  is  not  known  to  which  level  the  C02  (00^1)  level  is  Quenched,  although  the 
data  of  Cannemijer(24)  indicate  that  it  is  the  C02  (03*0)  level  in  C02-C02 
collisions  and  the  C02  (02f0)  level  in  C02-rare  gas  collisions.  However, 
even  if  it  is  the  near-resonant  04*0  or  other  nearby  levels  in  N20,  detailed 
balance  indicates  that  the  reverse  rate,  coupling  modes  vj  and  v  2  to  mode 
in  collisions  with  N2,  will  be  slow.  This  is  supported  by  the  experiments 
of  Yardley,  (2^)  who  excites  N20  mode  V3  in  mixtures  with  N2.  He  observes 
that  N20  mode  V3  rapidly  couples  to  N2  by  V-V  transfer  in  about  1100 
collisions.  However,  at  the  end  of  this  (short)  time  the  extent  of  residual 
excitation  in  mode  V3  indicates  that  modes  v]  and  v2  have  not  been  excited 
by  either  direct  intramolecular  V-V  processes  or  intermolecular  V-V 
transfer  back  from  N2.  (25)  Quenching  of  C02*  (v3)  to  modes  vi  and  y2 
by  H20  is  much  faster  than  by  N2,  requiring  on  the  order  of  103  collisions. 
This  probably  because  H20  rotation  can  absorb  much  of  the  nonresonance 
between  C02  (00°  1 )  and  C02  (mnf0).  However,  even  this  rate  of  coupling 
is  likely  to  be  much  slower  than  direct  quenching  of  modes  vi  and  v2  by 
H20  (see  below). 

We  are  therefore  led  to  neglect  the  coupling  of  modes  vj  and  v2  to 
mode  V3  and  to  assume  that  modes  vi  and  v2  are  rapidly  coupled  This 
leads  to  the  model  shown  in  Fig.  5.  According  to  this  model,  modes  ej  and 
v2  are  reduced  to  four  lumped  levels,  each  consisting  of  several  near¬ 
resonant  levels  as  indicated  by  the  dotted  lines. 
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Few  of  the  rates  of  relaxation  of  these  levels  have  been  measured 
experimentally.  The  relaxation  of  the  N20(vt)  model  has  been  determined 
by  acoustical  and  spectrophonic  techniques ( l ”)  and  found  to  be  1.8  x  10“^ 
cm^/s.  This  is  approximately  a  factor  of  five  faster  than  the  corresponding 
relaxation  rate  for  COolv^)’  as  roight  be  expected  from  the  slightly  lower 
vibrational  frequency  of  To  scale  the  rate  to  different  temperatures 

we  assume  that  the  temperature  dependence  for  ^0(^2)  +  N 2  is  the  same  as 
that  for  002(^2)  +  N2-  To  scale  the  rate  to  higher  levels  we  assume  the 
harmonic  oscillator  rule,  which  says  that  for  level  (0  m'O)  the  quenching 
rate  is  proportional  to  m. 

Continuing  the  analogy  to  CO2,  we  anticipate  that  the  quenching  of 
N20(v2)  by  H2O  should  be  a  very  rapid  process  since  the  quenching  of 
CC>2(v2)  requires  only  about  20  collisions.  (1?)  Two  theories  have  been 
advanced  to  explain  the  rapid  quenching  of  CO2.  Widom  and  Bauer(2b) 
suggest  that  CO2  and  H2O  interact  chemically,  since  H2CO3  is  known  to 
exist  in  solution.  Semiquantitative  calculations  based  on  this  mechanism 
give  the  right  order  of  magnitude  for  the  quenching  rate.  Sharma(^7) 
suggests  that  the  multipolar  interaction  between  CO2  and  H2O  is  responsible. 
Detailed  computations  based  on  this  mechanism  are  also  in  reasonable 
agreement  with  the  experimental  data.  If  the  latter  explanation  is  correct, 
a  similar  mechanism  should  operate  in  N2O-H2O  collisions.  On  this  basis 
we  assume  that  the  rate  of  quenching  of  N2O  (01  ^ 0 )  by  H2O  is  the  same  as 
the  measured  rate  of  quenching  of  CO2  f 0 1  ^ 0 )  by  H2O  at  all  temperatures. 

To  scale  to  higher  levels,  the  harmonic  oscilllot.'*-  rule  is  used. 

The  quenching  of  ^Ofvj)  is  likely  to  be  much  slower  because  of 
the  larger  '  ibrational  level  spacing,  compared  with  mode  v  2-  Provided 
that  modes  v  j  and  v-2  are  closely  coupled  the  direct  quenching  of  mode  v  1 
can  be  ignored.  However,  if  they  are  not  so  coupled  it  is  necessary  to 
consider  what  might  be  the  most  important  processes  relaxing  mode  vj. 

A  strong  candidate  for  the  dominant  role  is  near-resonant  V-V  transfer 
to  H2O. 

NzO  (10°0)  +  HzO  (000)  -*  NzO  (00°0)  +  HzO  (010)  -  310  cm"1  (21) 

Although  this  reaction  might  seem  far  from  resonant,  H2O  has  shown  itself 
capable  of  a  variety  of  not-so-resonant  processes,  apparently  because  a 
large  amount  of  nonresonance  can  be  made  up  by  the  H2O  rotational  fre¬ 
quency.  Similar  reactions  to  which  (21)  may  be  compared  include  the 
reactions 

CS  (v  =  1)  +  HzO  (000)  -  CS  (v  =  0)  +  HzO  (010)  -  323  cm'1,  (22) 

NO  (v  —  1)4"  HzO  (000)  -»  NO  (v  =  0)  +  HzO  (010)  +  281  cm"1,  (23) 

N20  (00°1)  t  HBr  (v  =  0)  -  NzO  (00°0)  +  HBr  (v  =  1)  -  335  cm'1  (24) 
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These  reaction  have  room  temperature  rates  (in  the  exothermic  direction) 
equal  to  1  x  10'11,  2x  10’12,  and  2  x  1 0“  1 2  cm3/s,  respectively.  ( 1  •)  A 
rate  equal  to  10-^2  cm^/sec  seems  reasonable  for  reaction  (21)  in  the 
exothermic  direction.  This  is  much  slower  then  the  quenching  of  mode 
v 2  by  H2O  and  will  be  negligible  if  the  modes  are  coupled.  However,  if 
modes  v  ^  and  are  not  coupled,  their  rates  of  relaxation  are  likely  to  be 
quite  different. 


These  processes  were  summarized  in  Fig.  5,  where  the  times 
indicated  correspond  to  the  atmospheric  conditions  of  Figs.  2,  3,  and  4. 
Only  a  small  amount  of  heating  accompanies  the  absorption  of  a  laser 
photon.  The  bulk  of  the  heating  occurs  on  a  time  scale  of  the  order  of 
2  x  10"^s,  corresponding  to  the  quenching  of  N^O  iv 2)  by  H^O. 


5.  Kinetics  ol  Absorption  of  DF  Laser  Radiation  by  CH4 

DF  laser  radiation  is  absorbed  by  CH4  on  the  V4  overtone  band 


CH4  (0000)  +  hv  -  CH4  (0002)  +  AE,  (25) 


and  the  combination  band 


CH4  (0000)  +  hv  -  CH4  (0101)  +  AE.  (26) 

-34) 

)  and 

less  than  100  collisions  to  couple  modes  and  V4).  (32)  since  this  coupling 
presumably  takes  place  via  intermolecular  V-V  transfer  in  high  concentra¬ 
tions  of  CH4,  the  modes  are  likely  to  be  much  more  slowly  coupled  in 
dilute  mixtures  of  CH4.  On  the  basis  of  his  experimental  data  Yardley 
concludes(32)  that  the  rate  of  coupling  by  rare  gases  is  probably  at  least 
a  factor  of  ten  slower.  Although  atmospheric  species  such  as  N2  and  O2 
probably  behave  more  like  rare  gases  than  like  CH4,  we  shall  nevertheless, 
assume  that  all  the  modes  are  strongly  coupled,  since  there  is  not  enough 
data  on  vibrational  relaxation  in  CH4  to  justify  a  more  detailed  description. 
Furthermore,  the  symmetric  and  asymmetric  bending  modes,  V2  and  V4, 
have  comparable  vibrational  frequencies  and  are  likely  to  behave  similarly, 
whether  coupled  or  not.  The  stretching  modes,  vi  and  V3,  have  much 
larger  vibrational  frequencies  and  probably  relax  much  more  slowly  than 
the  bending  modes.  They  will  therefore  have  little  effect  on  the  relaxation 
of  CH4  except  to  increase  the  effective  degeneracy  of  the  bending  levels 
to  which  they  are  coupled.  On  the  basis  of  these  considerations  we  adopt 
the  model  shown  in  Fig.  5. 


In  the  experiments  which  have  been  done  on  CH4  to  date  (18,  28 
all  the  vibrational  modes  have  been  coupled  very  rapidly  by  CH4-CH4 
collisions  (on  the  order  of  70  collisions  to  couple  modes  vt.  and  v  a. 
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The  relaxation  of  CH4  in  the  atmosj  ere  is  dominated  by  V-V  trans¬ 
its  to  O2  and  H2O.  The  rate  of  transfer  from  CH4  to  O 2  is  well  document¬ 
ed,  (31)  although  it  is  not  known  whether  it  lakes  place  from  near  -  resonant 
(0100)  level  or  the  infrared  active  (0001)  level.  Under  the  assumption  that 
these  levels  are  in  equilibrium  with  each  other  we  assign  this  process  a 
room  temperature  rate  equal  to  6  x  10"  13  cm^/s  in  the  exothermic  direction 
(O2  transfer  to  CH4).  The  reverse  rate  is  obtained  by  detailed  balance, 
using  the  partition  function  of  the  combined  CH4  (0001)  and  CH4  (0101)  levels. 
To  scale  this  rate  to  higher  levels  we  use  the  harmonic  oscillator  rule,  as 
before.  To  scale  this  rate  to  different  temperatures  we  assume  that  the 
reaction  is  sufficiently  near  resonance  for  the  probability  per  collision  to 
vary  as  T*l.  The  rate  coefficient  in  this  case  varies  as  T“  1/2. 

The  rate  of  V-V  transfer  to  H20(v>2)  is  not  so  well  documented,  and 
there  is  some  controversy  regarding  the  extraction  of  rate  coefficients  from 
the  limited  data  on  CH4/H2O  mixtures.  (1^>  30)  jf  we  accept  Bauer's  inter¬ 
pretation  of  the  data,  Tie  rate  of  transfer  from  H2O  (010)  to  Cfbifv2  +  V4) 
is  6  x  10"!  1  ern^/s  at  room  temperature.  To  scale  this  rate  to  higher 
levels  the  harmonic  oscillator  rule  is  used,  as  before.  To  scale  the  rate 
to  other  temperatures  we  assume  that  the  reaction  is  sufficiently  near 
resonance  that  the  rate  varies  as  T“l/2,  as  was  assumed  for  transfer 
to  O2.  The  subsequent  relaxation  of  H20(v2>  and  O 2  was  discussed  in 
Section  II-  B.  3. 

These  processes  are  summarized  in  Fig.  6,  where  the  relaxation 
times  indicated  correspond  to  the  same  atmospheric  conditions  as  in 
Figs.  2  to  5.  Only  a  very  small  amount  of  heating  or  cooling  occurs  at  the 
absorption  of  a  laser  photon  by  CH4.  At  sufficiently  low  H2O  concentra¬ 
tions  («  0.  2%)  some  cooling  will  occur  on  the  fast  time  scale  (~4  x  10-^s 
in  Fig.  6)  corresponding  to  transfer  to  O2.  Such  small  concentrations  of 
H2O  can  occur  at  high  altitudes,  since  H2O  is  lot  uniformly  distributed 
throughout  the  atmosphere.  Although  V-V  trai  fer  from  CH4  (V2  and  V4) 
to  H2O  (v  2)  3®  endothermic,  this  process  will  not  lead  to  cooling  because 
'he  H2O  (i*  2 )  vibration  is  quenched  on  the  same  time  scale  with  a  large 
release  of  heat.  Thus,  under  conditions  of  sufficiently  high  humidity 
(»  0.  2%  H2O),  which  ordinarily  prevail  at  low  altitudes,  the  laser  energy 
absorbed  by  CH4  will  lead  to  heating  on  a  fast  time  scale  (~7  x  10-^s) 
corresponding  to  transfer  to  H20  (v2)  and  quenching  of  the  I^O  (v2) 
(principally  by  H2O  itself).  The  energy  transferred  to  O2  will  appear  as 
heat  on  a  longer  time  scale  (~10"^s)  corresponding  to  the  relaxation  of 
O2  via  H2O  (v2)- 

6.  Kinetics  of  Absorption  of  DF  Laser  Radiation  by  HDO 

HDO  absorbs  DF  laser  radiation  on  the  fundamental  transition 


HDO  (000)  f  hv  —  HDO  (100)  +AE. 


(27) 
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Fig.  6  Kinetics  of  Absorption  of  OF  Laser  Radiation  by  CH4 


No  data  exist  regarding  the  relaxation  of  HDO  (v  j )  (or  v2  or  v»3, 
for  that  matter).  However,  preliminary  data  have  been  obtained  at  AERL 
regarding  the  relaxation  of  H20  (v  i )  in  collisions  with  H20.  >3^)  These 
data  indicate  that  modes  v  j  and  v2  (and  1/3  as  well)  are  coupled  at  every 
H20  -  H20  collision,  despite  the  fact  that  the  levels  are  separated  by  about 
500  cm-  1  Because  of  the  similarity  of  HDO  and  H20,  we  assume  that 
HDO  (vi)  and  (v2)  are  coupled  to  each  other  and  to  H20  (v2)  at  every 
collision  by  the  reactions 


HDC  (100)  +  HzO  -  HDO  (020)  +  H20  (28) 

H20  (010)  +  HDO  (000)  -*  HzO  (000)  +  HDO  (010).  (29) 

These  reactions  are  assumed  to  have  rate  constants  equal  to  5  x  10"  11 
cm3/sec.  To  scale  to  higher  levels,  the  harmonic  oscillator  rules  are 
used.  The  T-V  relaxation  of  HDO  (v2)  by  H20,  N2  and  02  is  assumed  to 
be  the  same  as  for  H20  (v2). 

In  addition  to  these  rapid  H20  -  HDO  processes,  near  resonant  V-V 
transfer  to  N2  and  02  may  be  important,  especially  under  conditions  of  low 
humidity: 

Oz  (v  =  1)  +  HDO  (000)  -  02  (v  =  0)  +  HDO  (010)  +  153  cm"1,  (30) 

HDO  (100)  +  N2  (v  =  0)  -  HDO  (000)  +  N2  (v  =  1)  +  393  cm'1.  (31) 

The  rate  of  reaction  (30)  is  estimated  to  be  about  1. 0  x  lO'1^  cm3/s.  This 
is  slightly  slower  (30%)  than  the  corresponding  rate  of  transfer  between 
H20  (010)  and  02,  which  is  more  nearly  resonant.  It  is  more  uifficult  to 
estimate  the  rate  of  reaction  (31).  Due  to  the  larger  nonresonance,  compared 
with  reaction  (30),  the  rate  of  reaction  (31)  has  somewhat  arbitrarily  been 
assumed  to  be  2  x  10”  * 3  cm3/s,  which  is  a  factor  oi  five  slower  than 
reaction  (30). 

These  processes  are  summarized  in  Fig.  7,  where  the  times  indica¬ 
ted  correspond  to  the  same  atmospheric  conditions  as  in  Figs.  2  to  6.  The 
dominant  effect  of  absorption  by  HDO  (vj  is  seen  to  be  rapid  heating  on  a 
time  scale  of  the  order  of  4  x  10~°s,  corresponding  to  the  quenching  of 
HDO  (v 2)  and  H20  (v2)  by  H20. 

C.  RELAXATION  EQUATIONS 

V/ith  these  approximations  it  is  a  simple  matter  to  write  down  the 
kinetic  equations  describing  the  relaxation  of  the  system.  For  this  purpose 
we  shall  refer  to  a  group  of  closely  coupled  vibrational  and  rotational  states 
as  "level  i",  and  define  the  partition  function  Qj  and  mean  energy  for 
level  i  by  the  expressions 
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Fig.  7  Kinetics  of  Absorption  of  DF  Laser  Radiation  by  HDO 
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where  gm  is  the  degeneracy  of  state  m,  Em  is  the  energy  of  state  m, 
k  is  Boltzmann's  constant  and  T  is  the  temperature. 


The  density  n^  of  molecules  in  level  i  satisfies  the  inviscid  equation 


(34) 


where  D/Dt  is  the  usual  convective  derivative,  and  v  is  the  local  velocity. 
Thus,  each  level  of  each  molecular  species  is  treated  as  a  chemical  species 
capable  of  undergoing  the  reactions  represented  by  the  right-hand  side  of 
the  equation.  The  first  term  on  the  right-hand  side  represents  the  effect 
of  vibrational  relaxation  processes  of  the  type 


A.  +  B.  -*  A.  +  B„, 
i  k  )  t 


(35) 


in  which  molecule  A  in  level  i  reacts  with  molecule  B  in  level  k  to  form 
molecule  A  in  level  j  and  molecule  B  in  level  £,  with  rate  coefficient 
k-*-  ( 

k  •  To  express  this  term  in  the  form  shown,  the  usual  detailed 
balance  relation 


QjQC 


e-k 


Q.  Q, 

i  k 
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has  been  used.  The  second  term  on  the  right-hand  side  represents  the 
effect  of  absorption  processes  of  the  type 


A.  +  hv  -  A.  (37) 
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and  stimulated  emission  processes  of  the  type 


A.  +  hv  A.  +  2hv.  (38) 

J  i 


In  Eq.  (34)  ct/VJ  is  the  cross  section  at  the  frequency  v  and  <5^  is  the  flux 

of  photons  at  this  frequency.  To  express  this  term  in  the  form  shown,  the 
detailed  balance  relation 


a 


(v) 


Qi  exp 


(-  — ) 
'  kT ' 


(39) 


has  been  used.  This  relation  is  derived  in  Appendix  A. 

To  solve  completely  the  problem  of  laser  beam  propagation  it  is 
necessary  to  solve  the  vibrational  relaxation  equations  together  with  the 
hydrodynamical  momentum  and  energy  equations  describing  the  motion  of 
the  gas  and  Maxwell" s  equations  (or  an  eiconal  approximation)  describing 
the  propagation  of  the  beam.  However,  to  simplify  the  study  of  absorption 
kinetics,  the  effects  of  beam  propagation  are  ignored,  and  the  laser  flux 
is  specified  as  a  function  of  time.  To  further  streamline  the  problem, 
the  dydrodynamic  effects  are  approximated  by  assuming  either  constant 
pressure  p  or  constant  total  density  n,  where 

-  ■£  ».• 

i 

p  =  nkT. 


(40) 

(41) 


The  constant  pressure  approximation  is  appropriate  when  the  laaer 
pulse  is  so  long  that  a  pressure  pulse  generated  by  the  heating  of  the  gas 
can  propagate  beyond  the  width  of  the  beam.  Thus,  if  W  is  a  characteristic 
width  of  the  laser  beam  and  At  is  the  pulse  duration  (as  seen  by  a  "particle" 
of  gas,  in  the  case  of  a  moving  beam)  then  the  constant  pressure  approxima¬ 
tion  is  valid  when 


aAt/W  »  1, 


(42) 


where  a  is  the  speed  of  uound.  For  a  sound  speed  a  3x10  cm/s,  and  a 
beam  width  W~10cm,  we  see  thst  the  constant  pressure  approximation  is 
appropriate  for  pulse  lengths  At  >  1ms.  The  constant  density  approximation 
is  appropriate  when  the  laser  pulse  is  so  short  that  a  pressure  pulse  has 
no  time  to  propagate  a  significant  distance,  so  that  the  gas  cannot  relax 
hydrodynamic  ally .  This  limit  is  valid  when 
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a  At/W  «  1. 


(43) 


Under  the  above  assumptions,  this  corresponds  to  pulse  longthsAt<  100  jis. 

In  the  constant  pressure  limit  the  energy  equation  may  be  expressed 

"  ~§T  (f  kT  +  €  )  hv  4>v  n.  o.  (v).  (44) 

v,  i 

where 

n?  ni  €i  (45) 


is  the  total  internal  (rotational  +  vibrational)  energy.  This  equation  must 
be  solved  simultaneously  with  the  relaxation  equations.  In  the  constant 
density  limit  the  energy  equation  may  be  expressed 

n  (f  kT  +  e)  hv  lv  n.  q.  (v).  (46) 

v,  i 

In  the  remainder  of  this  report,  only  constant  pressure  solutions 
will  be  discussed.  However,  the  temperature  change  AT  is  invariably 
small  (AT/T  «  1),  so  that  the  relaxation  equations  are  essentially  decoupled 
from  the  energy  equation.  Under  these  conditions  the  constant  density 
temperature  change  ATn  is  simply  related  to  the  constant  pressure  tempera¬ 
ture  change  ATp  by  the  formula 

AT  =  t  AT  .  (47) 

n  3  p  '  ' 


III.  RESULTS 


To  completely  assess  the  importance  of  kinetics  for  the  propa¬ 
gation  of  DF  and  CO  laser  radiation  through  the  atmosphere,  it  is 
necessary  to  survey  a  wide  range  of  laser  parameters  (especially  intensity 
and  pulse  length)  and  atmospheric  conditions  (especially  altitude  and 
humidity).  Although  this  task  has  not  been  completed,  preliminary 
results  have  been  obtained  and  these  are  discussed  below 

Figures  8  through  15  (interspersed  in  the  text  as  discussed)  give 
results  illustrating  the  use  of  our  analytical  model  for  atmospheric 
absorption  of  DF.  Results  at  CO  laser  wavelengths  are  given  in  Fig.  16. 
For  all  cases  shown  here  we  have,  for  simplicity,  used  a  square  wave, 
30-p.s  pulse,  and  assumed  the  same  power  level  on  each  laser  line  with 
all  lasing  occurring  at  the  same  time.  It  should  be  noted  that  the  model 
is  in  no  way  restricted  to  these  conditions. 

The  concentration  of  atmospheric  absorbers  is  given  in  Table  V. 
With  the  exception  of  water  vapor,  which  has  a  variable  concentration, 
the  concentrations  are  determined  by  the  natural  atmospheric  abundance. 
For  water  vapor,  we  have  used  two  limiting  concentrations  at  sea  level 
which  were  chosen  to  represent  very  dry  (2  torr  H2O)  and  very  wet 
(20  torr  H2O)  atmospheres  and  one  representative  concentration  at 
12  km  altitude,  5.  6  x  10~  3  torr  H2O.  (All  calculations  were  performed 
for  the  sea  level  temperature  of  300°K.  The  error  thus  introduced  at 
the  higher  altitude  is  probably  small.  As  stated  above,  absorption  by 
both  the  N2  and  H2O  continua  are  neglected.  ) 

Calculations  were  performed  at  two  altitudes,  sea  level  and 
12  km,  with  12  km  chosen  as  a  reasonable  maximum  altitude  of  interest. 


TABLE  V 

CONCENTRATION  OF  ABSORBING  MOLECULES 


Absorbing  Molecule 

Concentration 

N^O 

0.  28  ppm 

ch4 

1 .  6  ppm 

C02 

330  ppm 

H  *0  -  sea  level 

20  torr 

2  torr 

12  km 

5.  6  x  1 0~  ^  torr 

HDO 

3x  10-4  of  H2O  concentration 
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A.  ABSORPTION  OF  DF  LASER  RADIATION 


l.  Calculations  at  Sea  Level 

Figure  8  presents  the  translational  temperature  rise  as  a  function 
of  time  during  the  pulse  for  each  of  the  atmospheric  absorbers  separately 
anf  for  their  sum.  (It  should  be  noted  that  for  all  the  cases  we  have  run, 
the  temperature  rise  due  to  a  simultaneous  absorption  by  all  absorbers 
is  equal  to  the  sum  of  the  individual  temperature  rises.  )  1  he  results  of 

this  figure  represent  a  reasonable  upper  bound  in  that  the  incident  power 
was  taken  as  107  W/cm2  on  each  of  the  eight  laser  lines  and  the  water 
concentration  was  high  (20  torr).  (The  effect  of  water  concentration  on 
the  temperature  rise  is  discussed  below  and  may  be  obtained  by  comparing 
the  results  of  Fig.  8  with  those  of  Fig.  9.  )  Thermal  effects  after  the 
power  is  shut  off  are  not  shown  in  Fig.  8.  After  the  power  is  shut  off, 
the  temperature  continues  to  rise  at  a  reduced  rate  due  to  th  ■  continuing 
collisional  deactivation  of  stored  vibrational  energy.  The  major  absorber 
under  the  conditions  of  Fig.  8  is  HDO-H2O  and  the  effects  of  CO2  and 
CH4  on  the  temperature  rise  are  insignificant.  Although  it  is  not  obvious 
in  the  figure,  N2O  is  slightly  bleached  under  these  conditions;  i.  e.  ,  the 
temperature  rise  due  to  N2O  absorption  is  less  than  linearly  related  to 
the  incident  power.  For  the  other  absorbers,  which  are  not  bleached, 
linear  extrapolation  to  results  at  lower  power  levels  may  be  performed. 

Figure  9  presents  the  same  information  as  that  shown  in  Fig.  8 
for  conditions  in  which  the  H2O  concentration  is  2  torr,  a  reasonable 
lower  limit  approximating  very  dry  conditions.  An  obvious  difference 
between  Figs.  8  and  9  is  that  at  the  lower  H2O  concentration  in  Fig.  9 
the  total  temperature  rise  is  approximately  one  order  of  magnitude  less 
than  the  temperature  rise  under  the  20  torr  H2O  conditions  of  Fig.  8. 

The  effect  of  the  H2O  concentration  on  the  temperature  rise  is  twofold, 
since  H2O  is  both  a  direct  absorber  of  laser  energy  and  a  major  participant 
in  the  kinetics  of  T-V  deactivation. 

Because  H 2 ^  deactivates  both  excited  N2O  and  H2O-HDO  by  T-V 
processes,  a  higher  H2O  concentration  results  in  more  absorbed  energy 
appearing  as  translational  energy  during  the  pulse,  both  for  H2O-HDO 
and  N2O.  Comparing  Figs.  8  and  9,  it  is  seen  that  N2O  has  become  the 
major  absorber  under  the  conditions  of  2  torr  H2O.  This  is  true  even 
though  N2O  is  more  bleached  under  these  conditions  than  at  higher  humidity 
If  collisions  with  water  were  the  only  deactivation  path  available  to  excited 
N2O,  then  the  temperature  rise  of  Fig.  9  due  to  N2O  wou’d  be  0.  1  that  of 
Fig.  8.  The  presence  of  other  deactivation  paths  makes  the  difference 
between  the  N2O  temperature  rise  of  Figs.  8  and  9  less  than  a  factor  of  10 
Comparison  of  the  temperature  rise  due  to  H2O-HDO  shows  the  low 
humidity  results  to  be  between  10  and  100  times  lower  than  the  high 
humidity  results.  Thus  N2O  becomes  the  principal  absorber  under  condi¬ 
tions  of  low  humidity.  Here  again  t1  c  effe't  is  due  to  both  lowering  the 
absorber  concentration  (from  20  to  2  Lorr)  and  lowering  the  deactivator 
concentration.  If  collision  with  water  were  the  only  deactivation  path 
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Fig.  9  Translational  Temperature  Rise  for  Atmospheric  Absorption 
at  DF  Laser  Wavelengths  at  Sea  Level 


available  to  excited  H2O-HDO,  then  the  H2O-HDO  temperature  rise  of 

Fig.  9  would  be  0.  01  that  of  Fig.  8.  If  there  were  no  deactivation  path 

via  water,  then  the  temperature  rise  for  H2O-HDO  of  Fig.  9  would  be 

0.  1  that  of  Fig.  8.  The  fact  that  there  are  other  deactivation  paths  available, 

makes  the  difference  between  the  H2O-HDO  temperature  rise  of  Figs. 

8  and  9  less  than  a  factor  of  100  and  greater  than  10.  The  differences 
between  the  effect  of  changing  the  HzO  concentration  on  absorption  by 
CH4  and  CO2  are  due  to  the  differences  between  the  T-V  deactivation 
effects  of  H2O  on  these  two  absorbers.  These  differences  are  not  detailed 
here  because  neither  is  an  important  absorber. 

Figure  10  presents  the  translational  temperature  rise  as  a  function 
of  time  during  the  pulse  for  each  of  the  atmospheric  absorbers  separately 
and  their  sum  for  an  incident  power  of  10^  W/cm  on  each  line  and  a 
water  concentration  of  20  torr.  (These  conditions  are  the  same  as  those 
of  Fig.  8  except  for  a  factor  of  10  lower  incident  power.  )  Becaure  N2O 
is  only  very  marginally  bleached  under  these  conditions,  the  temperature 
rise  due  to  N2O  is  closer  to  the  total  temperature  rise  than  in  F.g.  8. 

Also,  the  N2O  produced  temperature  rise  may  be  extrapolated  linearly 
from  the  results  of  Fig.  10  to  lower  power  levels.  The  temperature  rise 
due  to  all  other  absorbers  is  a  factor  of  10  less  than  that  of  Fig.  8. 

Figure  11  presents  the  translational  temperature  rise  at  20  ps 
due  to  absorption  by  NzO  only,  as  a  function  of  incident  power  fo:  each 
of  the  seven  laser  lines  which  are  absorbed  by  N2O.  (Since  3P  (8)  and 
2P  (10)  have  the  same  cross  section  for  absorption  by  N2O,  they  are 
presented  as  one  line  only  resulting  in  six  apparent  lines  in  the  figure.  ) 

As  is  expected,  the  magnitude  of  the  temperature  rise  is  proportional 
to  the  N2O  absorption  cross  section  with  the  largest  AT  for  the  laser  line(s) 
with  the  largest  N2O  absorption  cross  section.  (These  results  at  10  W/cm 
are  slightly  different  from  those  of  Fig.  8  because  here  each  line  was 
considered  separately.  For  the  case  in  which  all  lines  are  considered 
together  (Fig.  8),  the  total  AT  is  somewhat  less  than  the  sum  of  all  lines 
taken  separately  since  absorption  on  all  lines  depletes  a  common  lower 
state. ) 


Figure  12  shows,  for  the  laser  line  of  largest  N2O  absorption  cross 
section,  the  power  at  20  s  into  temperature,  transmission,  absorption, 
and  stimulated  emission,  all  as  a  function  of  incident  power  for  20-torr 
H2O  concentration.  At  all  power  levels  the  largest  fraction  of  incident 
power  is  transmitted  with  the  absorption  in  a  l -km  path  ~  2%  of  the  trans¬ 
mitted  power.  When  the  change  in  net  power  absorbed  becomes  less  than 
linear  with  increase  in  incident  power,  the  gas  is  beginning  to  become 
bleached,  i.  e.  ,  the  lower  state  is  beginning  to  become  depleted  In  Fig  12 
this  effect  becomes  noticeable  at  an  incident  power  level  of  ~  10  W/crn^ 
on  the  individual  laser  line.  When  the  change  in  net  power  absorbed  (power 
absorbed  minus  stimulated  emission)  with  increase  in  incident  power 
approaches  zero,  the  gas  is  fully  bleached  and  any  further  increase  in 
incident  power  results  in  no  further  increase  in  temperature.  Complete 
bleaching  would  occur  in  Fig.  12  at  an  incident  power  level  significantly 
higher  than  those  shown  here. 
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Fig.  10  Translational  Temperature  Rise  for  Atmospheric 
Absorption  at  DF  Laser  Wavelengths  at  Sea  Level 
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Figure  13  presents  the  translational  temper  iture  rise  due  to 
absorption  by  CH4  at  220  pprn  and  2  ppm  under  the  same  conditions  as 
Fig.  8.  (The  natural  concentration  of  methane  in  the  atmosphere  is 
~  2  ppm.)  The  temperature  rise  due  to  DF  laser  absorption  by  CH4, 
which  occurs  on  only  3  of  the  8  laser  lines  chosen  for  this  study,  is  linear 
with  CH4  concentration  as  expected.  If  one  arbitrarily  assumes  that  the 
temperature  rise  due  to  absorption  by  some  man-made  or  natural  pollutant 
of  interest  can  be  approximated  by  methane,  then  the  results  shown  in 
Fig.  1^  indicate  that  somewhere  around  a  concentration  of  100  ppm  the 
pollutant  will  become  the  principal  absorber  of  DF  laser  radiation  under 
the  conditions  of  Fig.  8.  Of  course,  there  is  no  sufficient  reason  other 
than  convenience  to  assume  absorption  and  kinetic  characteristics  for 
the  pollutant  similar  to  those  of  CH4,  and  these  results  should  be  treated 
as  of  heuristic  value  only.  For  any  real  case  of  interest,  knowledge  of 
the  absorption  cross  section  and  kinetic  scheme  would  permit  the  correct 
calculation  to  be  performed. 

2.  Calculations  at  12  Kilometer  Altitude 

Figures  14  and  15  present  the  results  of  calculations  performed 
at  an  altitude  of  12  km  as  contrasted  to  Figs.  8  through  13  which  show 
results  for  sea  level  conditions.  An  altitude  of  12  km  was  chosen  to 
represent  a  reasonable  highest  altitude  of  interest.  The  water  concen¬ 
tration  was  chosen  to  be  5.  6  x  10-3  torr,  which  is  representative  of 
conditions  at  12  km.  For  calculations  at  12  km,  the  cross  sections  for 
absorption  were  approximately  scaled  by  the  decrease  ;n  pressure  from 
l  atm.  This  is  a  reasonable  procedure  because  for  all  the  lines  of 
interest,  the  laser  line  is  far  enough  away  from  the  absorption  line  center 
so  that  the  cross  sections  scale  as  the  half-width  which  is  linear  with 
pressure  for  a  Lorentz  shape  line. 

Figure  14  jresents  the  12-km  translational  temperature  rise  as 
a  function  of  time  during  the  pulse  for  each  of  the  atmospheric  absorbers 
separately  and  their  sum.  The  incident  power  was  taken  as  5  x  10&  W/cm^ 
on  each  of  the  eight  lines.  Even  though  N2O  is  more  bleached  than  at  an 
incident  power  level  of  10^  W/cm2  on  each  line  at  sea  level  (compare  slope 
of  AT  for  N2O  here  with  that  of  Fig.  8  and  also  compare  Fig.  15  w'ith 
Fig.  11),  N2O  is  the  principal  absorber.  Because  the  H2O  concentration 
is  so  low,  HDO-H2O  become  insignificant  contributors  to  the  translational 
temperature  rise.  The  CO2  induced  temperature  rise  is  almost  comparable 
in  magnitude  to  that  of  N2O.  The  effect  of  CH4  absorption  under  these 
conditions  is  to  cool  the  atmosphere  causing  a  temperature  decrease. 

In  all  cases,  the  initial  slope  of  AT  for  CH4  is  negative  because  V-V  transfer 
to  O2  and  H2O  is  endothermic  resulting  in  net  cooling.  Under  sea  level 
conditions  where  collisional  processes  cause  rapid  equilibration  this  is 
a  short-lived  effect  compared  to  pulse  length  and  net  translational  heating 
is  rapidly  achieved.  The  reduction  in  the  number  of  collisions  at  12  km 
results  in  the  cooling  effect  persisting  throughout  the  time  scale  of  the 
pulse.  As  is  always  true,  the  eventual  result  of  the  absorption  of  a  laser 
photon  is  heating,  but  at  12  km,  the  time  to  reach  equilibrium  is  longer 
than  the  pulse  time  for  CH^. 
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Fig.  13  Translational  Temperature  Rise  Due  to  Absorption  by 
CH4  at  DF  Laser  Wavelengths 
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Figure  15  is  the  12-km  analog  of  Fig.  11.  It  presents  the  transla¬ 
tional  temperature  rise  at  20  jisec  due  to  absorption  by  N2O  only  as  a 
function  of  incident  power.  The  increased  magnitude  of  bleaching  effects 
at  12  km  compared  with  the  sea  level  conditions  of  Fig.  1  1  is  easily  seen 
in  the  steep  decrease  in  slope  of  AT  as  higher  power  levels  are  achieved 

B.  ABSORPTION  OF  CO  LASER  RADIATION 


Figure  16  presents  the  sea  level  translational  temperature  rise 
as  a  function  of  time  during  the  pulse  for  absorption  by  H2O.  The  results 
shown  here  are  an  upper  bound  in  that  the  incident  power  is  taken  as  10^ 
W/cm^  on  each  of  the  12  CO  laser  lines  given  in  Table  III  and  the  I-LO 
concentration  is  20  torr.  Results  are  presented  both  for  the  AERL  best 
estimates  of  H2O  absorption  cross  sections  and  those  ofMcClatchey  ( 
The  AERL  results,  as  expected  from  Table  III,  are  almost  a  factor  of 
three  higher.  Since  there  are  no  bleaching  effects,  these  results  can 
linearly  extrapolated  to  lower  power  levels.  The  important  point  to  be 
noted  here  is  that  the  temperature  rise  due  to  absorption  at  CO  laser 
wavelengths  is  orders  of  magnitude  greater  than  at  DF  laser  wavelengths 
for  similar  conditions. 


TEMPERATURE  RISE  (°K) 


Fig.  16  Translational  Temperature  Rise  for  H2O  Absorption  at 
CO  Laser  Wavelengths  at  Sea  Level 
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APPENDIX  A 


DETAILED  BALANCE  FOR  OPTICAL  TRANSITIONS 


We  consider  a  molecule  N  whose  manifold  of  internal  levels  may  be 
grouped  into  two  sets  of  levels,  with  one  set  loosely  referred  to  as  the 
ground  state,  N  ,  and  the  other  loosely  referred  to  as  the  excited  state 
N  .  The  total  populations  of  these  states  are 


N°  =  £  N., 

(A-  la) 

n*  =  In., 
j  J 

(A-lb) 

where  N.  and  N.  are  the  populations  of  the  individual  levels  i  and  j, 
respectively.  We*  shall  assume  that  rapid  collisional  processes  keep  the 
relative  por'--1  ations  of  the  levels  within  each  state  in  equilibrium  with  each 
other  at  the  temperature  T.  That  is, 

g-exp(  -E./kT) 

N.  =  — - - - -  N°, 

1  Q° 

(A-Za) 

g  exp(-E  /kT) 

N.  =  _J - _J -  N  , 

J  Qr 

(A-2b) 

where  g.  and  g.  are  the  degeneracies  of  levels  i  and  j,  respectively 
and  E.  and  E.  4re  the  energies  (referred  to  the  same  zero)  of  levels  i 
and  j,1  respectively,  and  k  is  Boltzmann's  constant.  The  partition  func¬ 
tions  are  defined  by  the  expressions 

Q°(T)  =  £g.  exp(-E./kT),  (A-3a) 

i 

Q  (T)  =  £  g.  exp(-E./kT).  ( A-3b) 

j 

We  now  suppose  that  the  molecules  are  subjected  to  radiation  of  intensity 
Iv  and  frequency  v,  where 


hv  =  E.  -E., 
J  1 


(A-4) 
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and  h  is  Planck's  constant.  The  rate  of  induced  transitions  is  then 


N"  -N°  =  N.B  I  -  N.R..1,  (A-5) 

’  ij  j  ji 

where  Bjj  and  Bj^  are  Einstein's  coefficients  for  the  transitions  i-*j  and 
j"*i,  respectively.  These  coefficients  satisfy  the  relation 


g.B.. 

i  ij 


g.B... 
J  Ji 


(A-6) 


Substituting  Eqs.  (A-2),  (A-4)  and  (A-6)  into  Eq  .  (A-5)  we  obtain  the  re¬ 
sult 


N 


<.o  ,  _To  Q° 

V[  Q 


-hv/kT* 
e  '  N 


where  the  flux  of  photons  is 


(A-7) 


(A-8) 


and  the  absorption  cross  section  for  a  ground  state  molecule  is  defined 
by  the  expression 


CTvi<T>  = 


g^exp(-E^/kT 

Q° 


hvB. .. 
ij 


(A-9) 


It  shoqld  he  noted  that  if  several  transitions  are  absorbing  at  the  frequency 
v,  Eq.  (A-7)  may  be  used  simply  by  replacing  ov^  by  ov,  where 


CT 


v 


(A-10) 


If  several  frequencies  v  are  being  absorbed,  separate  equations  must  be 
written  for  each  frequency. 
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APPENDIX  B 


PARTITION  FUNCTIONS 


1.  NzO  "States" 

Referred  to  the  elemental  chemical  species,  the  partition  function 
(Q)  for  a  molecule  in  vibrational  level  ^^2^3)  is  given  by  the  expression 


V1V2  V3 


=  e 


xp(-AH  °AT)Q(T)Q>JVvR)^v 


(P-1) 


12 


where  AH  °  is  the  heat  of  formation  of  0°K,  k  is  Boltzmann's  constant 
and  T  is  me  temperature.  The  transitional  partition  function  {Q  J  is 
given  by  the  expression 

Q(T)  =  (27rmkT/h2)3/2,  <B‘2) 


where  m  is  the  molecular  mass  and  h  is  Planck's  constant.  According 
to  the  classical  rigid  rotor  approximation,  the  vibration-rotation  partition 
function  [Q<UR)]may  be  expressed 


Q(VR)C  =  g  kT/B  £ 
vlv2  v3  1  V1V2  V3 


X  exp 


[(Bvlv2£v 


r  -  g  £ 

3  V1V2  V3 


M 


(B-3) 


where  BVlV2£v3  is  the  rotational  constant  and  GVlv2£v3  is  the  vibra¬ 
tional  energy.  The  vibrational  degeneracy  is  given  by  the  formula 

g£  =  1  for  £  =  0, 

=  2  otherwise.  (B-4) 
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■ 


When  the  vibrational  le.els  are  grouped  into  closely  coupled 
"states",  we  obtain  the  following  partition  functions: 


Qo  Qoo°o 

(B-5a 

Q1  ""  Q01'0 

(B-5b 

Q2  =  Q02°0+Q0220  +  Q100 

( B  -  5c 

Q3  "  Q03'0  +  Q0330  f  Q1 1 '0 

(B-5d 

°4  Q04°0  +  Q0420  +  Q0440 

+  Q1 2°0  +  Q1220  +  Q20°0 

(B-5e 

2.  CH  "States" 

4 

Referred  to  the  elemental  chemical  species,  the  partition  function 
for  a  molecule  in  vibrational  level  (viv2v3v4^  is  givcn  the  expression 

°V1V2V3V4  =  EXp  Q(T)Q1^v3v4  (B-6) 

where  AHp°  is  the  heat  of  formation  of  0°K,  k  is  Boltzmann's  constant 
and  T  is  the  temperature.  The  translational  partition  function  is  given 
by  the  expression 


Q(T)  =  (2tt  mkT/h2)3/2 


(B-7) 


where  m  is  the  molecular  mass  and  h  is  Planck's  constant.  To  calculate 
the  vibration-rotation  partition  function  we  ignore  the  fine  multiplet  split¬ 
ting  of  the  intercombination  levels,  and  neglect  vibration-rotation  inter¬ 
action.  The  vibration-rotation  partition  function  may  then  be  expressed 


Q 


(VR) 

V1V2V3V4 


G 

v 


1V2V3V4 

kT 


) 
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(B-8) 


where  B  =  5.  24059  cm"*  is  the  rotational  constant  and  is  the 

vibrational  energy.  The  vibrational  degeneracy  is  given  by  the  formula 

gv  =  i  <v2+l)  (v4+l)  (v4+2)  (v3  +  1)  (v3+2),  (B-9) 

and  the  factor  12  accounts  for  the  rotational  symmetry. 

When  the  vibrational  levels  are  grouped  into  closely  coupled  "states", 
we  obtain  the  following  partition  functions: 


Qo 

=  Qoooo 

(B-lOa) 

Q, 

Q0001 

+  Q0100 

(B-10b) 

=  Q0002 

+  Q0101  +  Q0200 

(B-10c) 

+  Q1000 

+  Q0010 

(B-10d) 

ippipw-.-'p.-'-lT 


>  wipyufran  J,n  .,1.111 1*  u  (ill  I  UJI MWI  ji  n  wnyp^ 


APPENDIX  C 

SUMMARY  OF  VIBRATIONAL  RELAXATION  PROCESSES  AND  RATES 


Rate  k 

Process 

(cm  V sec) 

(1) 

N2(v=1)  +  HzO  —  N2(v=0)  +  H20 

-  1 5 

3.  6x10 

(2) 

N2(v=1)  +  N2  -  N2(v=0)  +  N2 

-  24 

2.  5x10 

(3) 

N2(v=1)  +  02  —  N2(v=0)  +  02 

2. 5xl0"24 

(4) 

N2(v=1)  +  O2(v=0)  —  N2(v=0)  +  02(v=l) 

1 .  6x10' 18 

(5) 

CO2(001)  +  N2(v=0)  —  C02(000)  +  N2(v=1) 

-  13 

5.4x10 

(6) 

co2(ooi)  +  h2o  —  co2(03o)  +  h2o 

-  1  3 

4. 3x10  ‘ 

(7) 

CO2(00l)  +  N2  —  C02(030)  +  N2 

-  1  5 

1.  9x10 

(8) 

CO2(001)  +  o2  -  CO2(030)  +  o2 

-  1  5 

1.  9x10 

(9) 

CO2(001)  +  C02  -*■  C02(030)  +  co2 

-  1  5 

4.  2x10 

(10) 

CO2(001)  +  HzO  —  C02(000)  +  h2o 

-  1  5 

4.  2x10 

(H) 

C02(030)  +  C02(000)  —  C02(020)  +  CO2(010) 

1. 0x10" 1 1 

(12) 

C02(020)  +  C02(000)  —  CO2(010)  +  CO2(010) 

6.  9x10"  12 

(13) 

CO2(100)  +  anything  -*  C02(020)  +  anything 

1.  Ox  1 0  1  ° 

(14) 

co2(oio)  +  h2o  —  co2(ooo)  +  h2o 

1.  8x1  0"  1  1 

(15) 

co2(020)  +  h2o  —  co2(oio)  +  h2o 

3.  6x10' 1 1 

(16) 

C02(030)  +  H20  -  C02(020)  +  HzO 

5.  4x1  0  1  1 

(17) 

CO2(010)  +  N2  —  C02(000)  +  N2 

3.  3x10"  15 

(18) 

C02(020)  +  N2  -*CO2(010)  +  N2 

7.  6xl0'15 

(19) 

C02(030)  +  N2  -~C02(020)  +  N2 

1.  1x10"  14 

(20) 

CO2(010)  +  02  —  C02(000)  +  o2 

3. 8xlO"1S 
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Rate  k 


Process 

(cmV  sec) 

(21) 

co2(020)  +  o2  ~co2(oio)  +  o2 

7.  6xl0'15 

(22) 

C02(030)  +  02  —  CO2(030)  +  02 

1.  lxlO"14 

(23) 

co2(oio)  +  co2  -  co2(ooo)  +  co2 

-  1  6 

7.  2x10 

(24) 

00,(020)  +  CO,  -  CO.  (010)  +  CO, 

-  1  4 

1. 4x10 

(25) 

C02(030)  +  C02  -  C02(020)  +  C02 

-  1 4 

2.2x10 

(26) 

H2O(010)  +  HzO  -H20  (000)  +  h2o 

4.  2xl0“n 

(27) 

H20(020)  +  HzO -H2O(010)  +h2o 

8.  3x10'  1  1 

(28) 

H2O(010)  +  N2  -«K20(000)  +  N2 

-  1  4 

6.0x10 

(29) 

H20(020)  +  N?  -H2O(010)  +  N2 

-  1  3 

1.  2x10 

(30) 

H2O(010)  +  02  -H20(000)  +  o2 

-  1  5 

2.  7x10  1 

(31) 

H20(020)  +  02  -H2O(010)  +  o2 

5.  4x10“  15 

(32) 

H2O(010)  +  O2(v=0)  -H20(000)  +  Oz(v=l) 

-  1  2 

1.5x10  1 

(33) 

H20(020)  +  O2(v=0)  -  H2O(010)  +  02(v=l) 

-  1  2 

3.  0x10  ' 

(34) 

02(v=l)  +  HzO  -  O2(v=0)  +  H?0 

3x1  0  *  4 

(35) 

02(v=l )  +  N2  —  O2(v=0)  +  N2 

3. 6xl0'19 

(36) 

02(v=l)  +  Oz  -  O2(v=0)  +  02 

3.  6x1  0_  1 9 

(37) 

o2(v=i)  +  co2(ooo)  - O2(v=0)  +  co2(oio) 

6.  6xl0"15 

(38) 

N2O(010)  +  H20  —  N20(000)  +  h2o 

2.  OxlO'1  1 

(39) 

N20(020)  +  H20  -  N2O(010)  +  h2o 

4. OxlO"1 1 

(40) 

N20(030)  +  H20-*  N2 0(020)  +  HzO 

6.  OxlO"1  1 

(41) 

N20(040)  +  H20  —  N20(030)  +  HzO 

8.  OxlO'1  1 

(42) 

N2O(010)  +  N2  -  N20(000)  +  N2 

1. 8x10"  ‘4 

(43) 

N2O(020)  +  N2  -  N2O(010)  +  N2 

-  1  4 

3.  6x10 

(44) 

N20(030)  +  N2  -  N2 0(020)  +  N2 

-  1 4 

5.  3x10  1 

-56- 


Rate  k 


Process 

(cmV  sec) 

(45) 

N20(040)  +  N2  —  N20(030)  +  N2 

-  1  4 

7.  1x10 

(46) 

n2o(oio)  +  o2  —  n2o(ooo)  +  o2 

-  1  4 

1.  8x10 

(47) 

N  ,0(020)  +  Oz  —  N2O(010)  +  Oz 

-  1  4 

3.  6x10 

(48) 

N20(030)  +  02  —  N20(020)  +  02 

-  1  4 

5.  3x10  1 

(49) 

N20(040)  +  02  —  N20(030)  +  02 

-  1 4 

7.  1x10 

(50) 

H2O(010)  +  N20(000)  —  H20(000)  +  N20(020) 

-  1  2 

1. 0x10 

(51) 

H2O(010)  +  N2O(010)  —  H20(000)  +  N20(030) 

1.  Oxl O' 12 

(52) 

H2O(010)  +  N2 0(020)  —  H20(000)  +  N20(040) 

2.  Oxl  O' 12 

(53) 

Oz(v=l)  +  CH4(0000)  —  O2(v=0)  +  CH4(0001) 

6.  OxlO-13 

(54) 

02(v=l)  +  CH4(0001)  —  O2(v=0)  +  CH4(0002) 

1.  2x1  O'  *  2 

(55) 

H2O(010)  +  CH4(0000)  —  H20(000)  +  CH4(0001) 

6.  4x1  0~  1  1 

(56) 

H2O(010)  +  CH4(0001)  —  H20(000)  +  CH4(0002) 

1.  3xlO'10 

(57) 

HDO(IOO)  +  HzO -HDO(020)  +  HzO 

2.  OxlO'10 

(58) 

HDO(OOO)  +  H2O(010)  — •  HDO(01  0)  +  H20(000) 

2.  OxlO'10 

(59) 

HDO(OIO)  +H20(010)  —  HDO(020)  +  H20(000) 

4.  OxlO'10 

(60) 

HDO(OOO)  +  H20(020)  —HDO(OIO)  +  H20(010) 

4.  OxlO"10 

(61) 

HDO(OIO)  +  H20  ^HDO(OOO)  +  HzO 

4.  2x1 0" 1  1 

(62) 

HDO(020)  +  H20  —  HDO(OIO)  +  H?0 

8.  3x10" 1  1 

(63) 

HDO(OIO)  +  N2  —  HDO(OOO)  +  N2 

-  1  4 

6.  OxlO 

(64) 

HDO(020)  +  N2  —  HDO(OIO)  +  N2 

-  1  3 

1.  2x10 

(65) 

HDO(OIO)  +  02  —  HDO(OOO)  +  02 

- 1  5 

2.  7x10 

(66) 

HDO(020)  +  02  —  HDO(OIO)  +  02 

- 1  5 

5.4x10 

(67) 

HDO(100;  +  N2(v=0)  —  HDO(OOO)  +  N2(v=l) 

- 1  3 

2.  Ox  i  0 

(68) 

02(v=l)  +  HDO(OOO)  —  O2(v=0)  +  HDO(OIO) 

1.  OxlO'12 

(69) 

Oz(v=l)  +  HDO(OIO)  —  02<v=0)  +  HDO(020) 

2.  OxlO 
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